A Possible Fermi liquid in the lightly doped Kitaev spin liquid 
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In this paper, we study the lightly doped Kitaev spin liquid (DKSL) and find it to be a Fermi 
liquid state. The DKSL has a closed Fermi surface around the center of the Brillouin zone with the 
quantized volume determined by the Luttinger's theorem. It has well-defined quasiparticles near 
the Fermi surface and can be detected in the angle-resolved photoemission spectroscopy. The DKSL 
has the topological Kitaev spin liquid (KSL) as its neutral background. It violates the Wiedemann- 
Franz law and has a large Wilson ration. These results has the potential experimental test in the 
iridates upon doping. 



Doping a spin liquid brings new phases in the strongly 
correlated electron systems and is believed to be the key 
physics in the high T c cuprates [JJ. Actually, the par- 
ent compounds of the cupratcs arc antifcrromagnetic or- 
dered, not in a spin liquid state. Many efforts have been 
made to find the spin liquid state in the quantum frus- 
tration Mott insulators[5]. In Ref. [31 Kitaev proposed an 
exactly solvable spin model with the spin liquid ground 
state on the honeycomb lattice. Recently the Kitaev 
model has the potential experimental realization in the 
strong spin orbit coupling magnets, such as the layered 
iridates A 2 Ir0 3 (A=Na, Li)[U [5]. In these materials, 
it is proposed in Refs. [5] and [7] that the spin interac- 
tions include both the isotropic Heisenberg term and the 
anisotropic Kitaev term. The Kitaev spin liquid (KSL) 
phase is stable when the Heisenberg term is small [3 [8]. 
The exact KSL state and its experimental realization pro- 
vide us the unique theoretical opportunity to test the 
physics of the doped spin liquid. This may give us some 
hints in the understanding of the high T c cuprates. 

The doped Kitaev spin liquid (DKSL) was studied in 
ReflHlby using the SU(2) slave boson method [Tl ITU] and 
in RefJTTlbv using the U(l) slave boson method (Tl?]. Both 
of them obtained the p-wave superconductors upon dop- 
ing. Motivated by these studies, we will make another 
investigation into the conducting state in the DKSL. In 
contrast to the p-wave superconductors in Refs 191 and ITT1 
we find a Fermi liquid state in this paper. Before dop- 
ing, the KSL has the topological order; it has four-fold 
degeneracy on the torus. The topological order is pro- 
tected by the Z 2 gauge structure and is robust against 
any local perturbations|13j- The KSL is a resonating 
valence bond (RVB) state [Til [T5] and the DKSL is rel- 
evant to high T c cuprates in the broader context of the 
doped RVB state[TJ [14]. In Refs. M and [TTJthe topo- 
logical robustness of the KSL breaks down after doping. 
None of the slave methods for the electron decomposi- 
tion satisfies the Z 2 gauge structure. To protect the Z 2 
gauge symmetry upon doping, we will use the dopon rep- 
resentation of the electron operators [TBI [17] to study the 
DKSL. The dopon theory is a full fcrmionic decomposi- 
tion of the physical electron operator. It describes the 




FIG. 1. The honeycomb lattice and the x, y and z links 
joining A and B sublattices. L = + |e2 and I2 = tb are 

the primitive vectors. 

DKSL in terms of two different fermionic components: 
"spinons", the neutral spin-1/2 excitations of the KSL; 
"dopons" , which describes the dopant holes and has the 
charge e and spin-1/2. In the lightly DKSL, the hybridiz- 
ing between spinons and dopons vanishes on the mean 
field level. The dopons form a Fermi liquid surrounded 
by the background spinons. The DKSL has well-defined 
low energy quasiparticles on the Fermi surface enclos- 
ing the quantized volume determined by the Luttinger's 
theorem [18]. The Fermi surface is electron like regard- 
less of whether we dope the holes or electrons into the 
KSL. These properties can be tested in the angle-resolved 
photoemission spectroscopy (ARPES) experiments. The 
background neutral spinons also have physical observ- 
able contributions. So unlike the standard Landau Fermi 
liquid [19] . the DKSL has a temperature-dependent spe- 
cific heat coefficient and violates the Wiedemann-Franz 
law in the transport measurements. In the T=0K limit, 
the DKSL has a large Wilson ratio. 

The KSL has the parent Hamiltonian of the spin S = 
1/2 system on the honeycomb lattice [3]: 

H K = -J ]T S? n S« 7 a = x,y,z (1) 

a-links 

where is a-component of the S = 1/2 spin operator 
on the site m and the summation runs over all the nearest 
neighbor a-links (a = x, y, z) joining A and B sublattices 
oriented in the a-th direction as shown in Fig. [TJ The 
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primitive vectors in the honeycomb lattice can be chosen 
as li = + \e 2 and l 2 = e 2 . 

To exactly solve the spin model ([T]), Kitaev introduced 
four Majorana fermions 6^, 6 TO , b z m and c m to rewrite 
the spin operators as = ib^Cm (a — x,y,z) with 
the constraint on every site D m = b x r fi v m b z m c m — \. In 
this paper, the Majorana fermions are normalized such 
as nrP\ — b t) a with sv x <y< z — b x - v ' z ^° — r Tn 

terms of the Majorana fermions, the Kitaev model ([I]) 
now has the form 



H 



K 



J 



a-links 



ic rn c n u rnn ^ ol y, z 



(2) 



with = ic£c£ and w m „ = -u* m . Without loss 

of generality, we assume conventionally that m is in A 
sublattice. Using the four Majorana fermions, we can 
also construct the fermionic Schwinger representation for 
the spin operators [H]: 



F„ 




a— x,y,2,0 



tv[(F m Fl-I)a%, 



x,y,z 



(3) 
(4) 



where a a (a — x,y,z) are the Pauli matrices and a° — 
il. We have the single-occupancy constraint for the 
Schwinger fermions G m = -|tr[( J F^F m - I)a a ] = 0. 
The mapping between the Majorana fermions and the 
Schwinger fermions is not unique and all such mappings 
are equivalent under SU(2) gauge transformations. The 
key observation of Kitaev is [H, w mn ] = and u mn are 
the constants of motion with eigenvalues U mn = ±1/2. 
The KSL is a RVB state with the parent mean field 
Hamiltonian [151 



Hk —J J]] U"„ii® 

a — links 



mn run 



2R 



vp 



U x u 



rnn ^mn 



x-links 



y- links 



z- links 



fllfll) + h.c, = K/^iW - + h.c and 

Km = l(/mf/nt + /rU4f) + h - c - Here ^vp is the energy 
of the nearest neighbor vortex-pair on the honeycomb 
lattice. [20] The Hamiltonian Q is invariant under local 
Z 2 transformations 



fn 



fmcr G m f^ 



u° 



where G m is an arbitrary function with only the two val- 
ues ±1. Such a Z 2 gauge symmetry leads to the topo- 
logical order with four-fold degeneracy for the KSL on 
the torus. This four-fold degeneracy is protected by 
the Z 2 gauge structure and is robust against any local 
perturbations [T5] . The wavefunciton of the KSL is a pro- 
jection of the ground state of the parent Hamiltonian ([5]), 



| KSL) = 7 7 |W)mf (P removes the double occupancy). 
Such a wavefunciton of the KSL is an exact result after 
the projection. 

Upon doping, we study the doped Kitaev model H = 
H t + Hk with the doping level x on every site. H t is the 
hopping Hamiltonian upon doping 



H t = -t J2 Vc \naCnaV 
(mn) a 



(7) 



We will employ the dopon representations for the electron 
operators [ini E] 
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T^dfrncr ( / . ^ frncr'd m — c 



(8) 



Here the spinon f ma and the dopon d ma are both 
fermionic operators. On every site, the states | f 0), 
| I 0), and the local singlet state |(| t|) — | 4-t)) are 
the physical states mapping onto the states | f), | J.) 
and the vacancy state |0), respectively. The operator 
Vd is to project out the unphysical triplet states be- 
tween the spinon and the dopon on every site. The 
spinon sector is always half-filled on the honeycomb lat- 
tice, ^2 a fmtr fma = L 

Here we slack our steps to give some explanations for 
such a dopon decompositions ^ . The dopon theory de- 
scribes the DKSL in terms of two fermions: "spinons" , 
the neutral spin- 1/2 excitations of the KSL; "dopons", 
which has the same charge e and spin-1/2 as the dopant 
holes. Here we start from the Mott insulator at half- 
filling and described the DKSL as a doped Mott insula- 
tor. 

The spinons in the dopon theory are always at half- 
filling on the honeycomb lattice. The Mottness is always 
described by the half-filled spinons even upon doping. 
The KSL at half-filling has the Z 2 gauge symmetry ^ 
and the nontrivial topological order. The topological or- 
der plays very important roles in spin liquids [T5] and ro- 
bust against local perturbations. However, the topolog- 
ical robustness is lost in the slave boson methods used 
in Refs. |9] and [H] even in the very lightly doped case. 
Because the physical electron operator Cj CT is no longer 
invariant under the Z 2 gauge transformation ([6]). This 
is very unacceptable if we have the simple belief in the 
topological order. Fortunately, the dopon theory still has 
the Z 2 gauge symmetry ([6| upon doping. In the dopon 
representation, the electron operator has two spinon /mo- 
multipliers and remains unchanged under the Z 2 gauge 
transformations. 

The dopons in the dopon theory are fermions differ- 
ent from the slave boson methods in Refs. l9landfTTl As 
noted in Ref|9l Willans et al showed in Ref. |2T|that a spin 
vacancy with a 7r-flux is stable as the low energy excita- 
tion in the dilute limit vacancy doped KSL. The doping 
process introduces the fermionic low energy excitations 
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for the charged dopant holes. So it is reasonable for us 
to treat the dopons as fermions in the lightly DKSL. 

The dopons have the spin- 1/2 degree of freedom which 
is another vital difference from the slave boson meth- 
ods in Refs. [§] and [TTJ The spinons are always at half- 
filling. So the spin degree of freedom in the dopon sector 
is used to neutralize the spin on the spin vacancy. Ref. 
I2T1 showed that the dopant holes have the low excitations 
as a spin vacancy with a 7r-flux. The flux in the KSL has 
something to do with the spin excitations described by 
in Eq. [ZU\ The spin- 1/2 degree of freedom in 
dopon sector is consistent with the numerical results. [2Tj 

Honestly, there is no rigorous proof for the validity of 
the dopon theory; nevertheless, it quite makes sense for 
us to use it in the lightly DKSL. We rewrite the hopping 
Hamiltonian as 

' V d [tT(a z D m Dl) - t r (D m a z (F m Fl - I) 



Hi 



(mn) 



xa z Dl) - tv(D m a z (F n F^ n - 1)g z D\) + tr(D m a z 
x(F m Fl - I)(F n Fl - l)o z D\a z )\V d (9) 



with the definition _D, 



d m 1; 



The process of 



hybridizing spinons and dopons, B m oc (tr(D m F m )) 7^ 0, 
breaks both the Z 2 gauge symmetry and the physical 
U(l) electromagnetic svmmetrv|161 IT?] resulting in the 
p-wave superconductors similar to the results in Refs. [5] 
and [TTj . However, the KSL has a gap E vp in the spec- 
trum. The spontaneous mixing B m competes with the 
gap energy E vp and is not allow for lightly doped case 
x < ~ 0.134. In this paper, we set t = J = 1 and 
E vp = 0.0668 J. 0] 

Without the mixing between the dopons and spinons, 
B m oc (tr(fl m f m )) = 0, the four-operator and six- 
operator terms in Eq. ([9| are decoupled into bilinear 
terms such as D m D n and F m F n on the mean field level. 
Therefore, the spinons and dopons are decoupled on the 
mean field level. The DKSL is described as the Fermi liq- 
uid formed by the dopons surrounded in the KSL formed 
by spinons. The full effective mean field Hamiltonian is 
now 



Hmf — Hd + Hf 
Here the dopon Hamiltonian reads out 



H d = 



£ ^k d I/k C r^k - + 2xN/J, a 



(10) 



(11) 



with£/ = 1,2 and ef 2k = ±§|e llt l1 +e lk ' 12 +1| -fj. d . N is 
the number of the unit cells. The chemical potential \x d 
is determined by the particle number of the dopant holes 
x. For x = 0.1, fj, d — —1.324. The spinon Hamiltonian is 

now(C/« n = i) 

Hf = J2'^U M ^kt + ^U m ^m (12) 




FIG. 2. [Color online] (a) The ARPES intensity J(k, 0) near 
the Fermi energy for the DKSL at a; = 0.1. (b) The disper- 
sions in the ARPES measurements. 



Here > takes the summation over half the Brillouin 
*-~> k 

zone. VL = (/Lr /La /l-k" > CT =t /I, is thc 

Nambu representation. The Hamiltonian matrices Mk CT 



/ z 5 _(k) -ig+(k)\ 

-ig-(-k) ig+(-k) 

-«5+(k) «5-(k) 

\ig+(-k) -t 5 _(-k) 

/ -ih+(k) -ih-(k)\ 

ih + (-k) ih-(-k) 

-ih-Qt) -ih + (k) 

\i/i_(-k) ih + (-k) 



(13) 



Here we have the definition g±(k) = ^(e 



ikl 2 



1) ± ^ and h ± (k) = ^-(e-*- 11 ± 
with the effective constants J e s = J[(l — x) 2 
and E^ p = E vp [(l — x) 2 — xt/2] . The spinon Hamilto- 
nian has the eigenvalues ±e{ k ^, ± e 2kf -^lkj. an( l = ' =e ik4. 
with the definition e{ kt = ia|e ikl1 -I- e lkl2 + 1| and 

fc 2kt ~~ fc lk^ — fc 2k^ — -°vp- 

In the dopon representation, we can write down the 
electron Green function G c in the DKSL 



-ik-li 



,-ik-l 5 



) 

ct/2] 



G c (k, w) oc G d (k, u) + Ginc 

= y 1 



(14) 



where G d (k, w) is the Green function for the dopon sector 
and Gi nc is the incoherent part. The dopon sector forms 
a Fermi liquid and the real part of the electron Green's 
function ReG c satisfies the Luttinger's theorem [T81 [T9"] 



i^ReG c (k,0) 

ktr 



> = 2x 



(15) 



ReG c (k, uj) is thc real part of the electron Green func- 
tion. The image part of the electron Green function is 
proportional to the intensity I(k, uj) oc ImG c (k, uj) in the 
ARPES measurements. The intensity at Fermi energy is 
shown in Fig. ga). The dispersion is shown in Fig. p2[ 
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FIG. 3. [Color online] The estimation of the temperature 
dependent Lorentz number for the DKSL at the doping x = 
0.1. 



Regardless of whether we dope the holes or electrons into 
the KSL, the Fermi surface is always electron- like pockets 
around the Brillouin center. 

The DKSL has well-defined low energy excitations and 
satisfies the Luttinger's theorem. It is a Fermi liquid. 
However, it has some unusual properties due to the exis- 
tence of the neutral spinon KSL in the background. The 
DKSL has the specific heat coefficient j — jd + lf, with 
the dopon specific heat coefficient 



7d 



= jE(^k) 2 /3sech 2 (/3^ k /2) 



(16) 



•c/ker 



and the spinon sector has the mean field specific heat 
coefficient 

V = 2 E'(^kJ 2 /?sech 2 (/3^ kff /2) (17) 

The specific heat coefficient in the DKSL is no longer 
temperature independent at low temperatures. In the 
transport measurements, the spinon sector has the con- 
tribution in the thermal conductivity, but no electric con- 
ductivity due to its charge neutrality. Therefore, the 
Wiedemann-Franz law breaks down. We can estimate 
the Lorentz number of the DKSL as follows 



K 

To- 



ld + 7/ 
Id 



(18) 



with L = ^ 



tj- (^r) ■ Here we assume that the 
mean free paths for the dopon and spinon are close to 
each other. More detail of the mean field path is beyond 
the scope of this paper and left to further investigation. 
For x — 0.1, the temperature dependent Lorentz number 
is shown in Fig. [3] There is a maximum at the temper- 
ature around T* = 0.31-E^. 

The spin susceptibility of the DKSL is given as x — 
Xd + Xj i with the dopon spin susceptibility 



(19) 



■E/ker 



which has the temperature independent Pauli behavior 
at low temperature and is proportional to the density of 



states in the dopon sector. The calculation of the spinon 
spin susceptibility is a little tricky [22]. The spinon 
sector is a p-wave superconductor and has vanished den- 
sity of states at zero temperature, however, the magnetic 
spin response is finite when the magnetic field is along 
the z direction. The mean field spin susceptibility in the 
spinon sector is given as 



with the definition 
1 

Okt = — } — 

e lkt £ 2kt 



«ki 



a a 



f f 

e lki e 2ki 



Xf 



fc lkt 



fc lk^ 



Off + a i 



f 4 g _(k) g _(-k) 1 

£ 2kt f J J 

fc lkt ' fc 2kt 

, 4Mk)M-k) 

£ 2k^ 



(20) 



fc lk^ 



/ 
2k4. 



k 



«ko 



(21) 



At finite doping x, the spinon spin susceptibility is given 

as 

J 



X/O) 



-x/(o) 



(22) 



For finite doping, the effective exchange constant J e g 
is reduced, however, the spinon magnetic response will 
be enhanced. This enhancement is also observed in the 
numerical results for the Kitaev model with the spin 
vacancy [23]. 

At zero temperature, the spinon sector has a vanished 
specific heat coefficient, but finite spin susceptibility. So 
it has a the Wilson ration 



R = 



Xd + Xf 



(23) 



Id + 7/ 

larger than 1 on the mean field level. For x = 0.1, R ~ 8. 

In this paper, we study the doped Kitaev spin liq- 
uid which has the potential experimental realization in 
the doped layered iridates. We calculate the electron 
Green's function and the thermodynamic properties for 
the DKSL which can be measured in the further exper- 
iments. We find the DKSL to be a Fermi liquid state 
but with a temperature dependent Lorentz number and 
a large Wilson ratio. 
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Supplementary material 



The spinon susceptibility 



In the following, we will calculate \f a t zero tempera- 
ture. To calculate the spinon magnetic spin response, we 
add the Zeeman term to the Kitaev bilinear model 



H f [h] =H S + \Y. ^L(M/ + AM h ^ 



(24) 



with Mh = diag(l, 1,-1,-1) and the chemical potential 
/if is to enforce single-occupancy constraint on the mean 
field level. At h — 0, [if = 0. The spinon sector now has 
the ^.-dependent eigenvalues ±e{ k ^(/i), ±e^ k ^(h), e{ k ^(/i) 
and el^ (h) with the definition 



"lkt 



-lkf. 



0) = \fpv 



7k, 



c 2kt 



0) = \J~Pv 



7k, 



(h) 



where 
Pk = 



^((4t) 2 + (4t) 2 ) + (M/ + ^) 2 , 



( e lkt) 2 ( e 2kt) 



7k = ( 
'k 



/ ^2 



-) 2 +4. 9 _(k) ff _(-k)( Ai/ + -) 2 , 



' " / \2 , (J \2\ , I,. _ h \2 



2 (Kk,r + (4k,) 2 ) + (M/-2^ 



(4 ) 2 - (4ki) 



*k = (- 



/ \2 



-) 2 + 4h + (k)h + (-k)([lf ~ - f 

The mean field free energy can be obtained from the 
mean field Hamiltonian 



Ff(h) 



tx{-pH s {h)) 



= -^E' ln (l + cosh[/3^ kCT (/ l )]) (25) 

P . 



■£/k<r 



The field dependent chemical potential is determined by 
dF f (h) _ 



= E^anM/?^)^) 9 ^ =0(26) 



At zero temperature, tanh(/3e) = 1. We take the partial 
derivative of Eq. ( [26] ) with respect to magnetic field h 
and then we have the following relation at h = 



E 



/ d 2 e f 



■e/k, 



with the definition 
1 

"kt 
«k4. 



ktr 



(27) 



/ / L c lkt ~ c 2kt 



e lkt £ 2kt 



J 



J 



f f l^lki ^ c 2k^ 
e lk^ e 2H 

Thus at h — 0, we have 

1 a t - a>i 



4 g -(k)g_(-k) 
/ , / ' 

e lkt + £ 2kt 

4fe + (k)/i+(-k) 

e lkl + e 2ki 



dh 



2«t + a i' 



(28) 



(29) 



The spin susceptibility is given as 




FIG. 4. The spin susceptibility at T = 0. 
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d 2 Fr ^-v / 1 9e{^i. t , 2 At /i = 0, we have a y, kg = and the spinon mean held 



g/j2 I 2 9ft ^"-h^j^W spin susceptibility at zero temperature 

.E/ker \ 

^ t ^ t . /2 „ ,30) ^-C^ + ^ + clf-^-^Pi) 



dh 2 



Also the spinon spin susceptibility Xf ls sensitive to the 
nearest neighbor vortex pair energy E vp . Without dop- 
ing, x = 0, the -E V p-dependent spinon spin susceptibility 
Xf is shown in Fig. [4] 



